How to determine an effective potential for a variable cosmological term 
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It is shown that if a variable cosmological term in the present Universe is described by a 
scalar field with minimal coupling to gravity and with some phenomenological self-interaction 
potential V((p), then this potential can be unambiguously determined from the following 
observational data: either from the behaviour of density perturbations in dustlike matter 
component as a function of redshift (given the Hubble constant additionally), or from the 
luminosity distance as a function of redshift (given the present density of dustlike matter in 
O | terms of the critical one) . 

<M ! PACS numbers: 98.80.Es, 98.80.Hw 

It has been known for many years that the flat Friedmann- Robert son- Walker (FRW) 
cosmological model with cold dark matter, a positive cosmological constant A > (fi +^A — 
1) and an approximately flat (or, Harrison-Zeldovich-like, n$ ~ 1) spectrum of primordial 
scalar (adiabatic) perturbations fits observational data better and has a larger admissible 
region of the parameters (Hq, Qq) than any other cosmological model with both inflationary 
and non-inflationary initial conditions (see, e.g., |2|). Here H is the Hubble constant, 
Qo = 87rGp m /3HQ includes baryons and (mainly) non-baryonic dark matter, = A/3Hq 
and the light velocity c = 1. This conclusion was based on the following arguments: a) 
relation between Hq and the age of the Universe to, b) the fact that observed mass/luminosity 
ratio never leads to values more than Q ~ 0.4 up to supercluster scales, c) comparison of 
cosmic microwave background (CMB) temperature anisotropies, power spectra of density 
. ^ ; and velocity matter perturbations, present abundance of galaxy clusters with predictions of 

^ ■ cosmological models with inflationary initial conditions; d) observed values of pt/ p m in rich 

galaxy clusters confronted with the range for the present baryon density pb admitted by the 
theory of primordial (Big Bang) nucleosynthesis. I don't include gravitational lensing tests 
(e.g., a number of lensed quasars) here, since conclusions based on them are less definite at 
present; however, the most recent reconsideration [[| has also led to a low value of Qq ~ 0.3. 

During last year two new pieces of strong evidence for Q < 1 have appeared. The 
first (historically) of them is based on the evolution of abundance of rich galaxy clusters 
with redshift z [§], see also the more recent paper || where the value Qq ~ 0.5 ± 0.2 (lcr 
uncertainty) is presented. Still, it should be noted that there have been already appeared 
some doubts on validity of the conclusion that Qq = 1 is really excluded ||. Much better 
observational data expected in near futute will help to resolve this dilemma unambiguously. 
The second, completely independent argument for Q = (0.2 — 0.4) follows from direct 
observations of supernovae (type la) explosions at high redshifts up to z ~ 1 [J/J. On the 
other hand, no direct evidence for a negative spatial curvature of the Universe (i.e., for the 
open FRW model) has been found. Just the opposite, the latest CMB constraints (based 
mainly on the results of the Saskatoon and CAT experiments) || , galaxy abundance at high 
redshifts [[J and the most recent analysis of the SNIa data in terms of an effective equation 
of state of a component adding Q to unity [jTOj strongly disfavor the open CDM model 
without a positive cosmological constant. Of course, the possibility to have both a positive 
cosmological constant and spatial curvature of any sign is not yet excluded, but, according 



to the "Okkam's razor" principle, it would be desirable not to introduce one more basic novel 
feature of the Universe (spatial curvature) without conclusive observational evidence. In any 
case, in spite of many theoretical and experimental attempts to exorcize it, a A-term is back 
again. 

It is clear that the introduction of a cosmological constant requires new and completely 
unknown physics in the region of ultra-low energies. Solutions with a cosmological constant 
occur in such fundamental theories as supergravity and M-theory. However, this cosmological 
constant is always negative and very large. As compared to such a basic "vacuum" state, a 
very small and positive cosmological constant allowed in the present Universe may be thought 
as corresponding to the energy density Ea of a highly excited (though still very symmetric) 
"background" state. So, it need not be very "fundamental". But then it is natural to omit 
the assumption that it should be exactly constant. In this case the name "a cosmological 
term" (or a A-term) is more relevant for it, so I shall use this one below. The principal 
difference between two kinds of non-baryonic dark matter - dustlike CDM and a A-term - is 
that the latter one is not gravitationally clustered up to scales ~ 30 hr 1 or more (otherwise 
we would return to the problem why Q observed from gravitational clustering is not equal 
to unity). Here h = H /100 km/c/Mpc. 

On the other hand, there exists a well-known strong argument showing that a A-term 
cannot change with time as fast as the matter density p m and the Ricci tensor (i.e., oc t~ 2 ) 
during the matter dominated stage (for redshifs z < 4 ■ 10 4 h 2 ). Really, if e\ oc p m , so that 
Q\ = const, then matter density perturbations in the CDM+baryon component grow as 
6 = (&) oc t a oc (1 + z)~ 3a/2 , a = V25-240 A -i ^ Ag a consequence; the total growth of 

perturbations A since the time of equality of matter and radiation energy densities up to 
the present moment is less than in the absence of the A-term. If Q\ <C 1, then A(Oa) = 
A(0)(1 — (6.4 + 2lnh)Q\). Since parameters of viable cosmological models are so tightly 
constrained that A may not be reduced by more than twice approximately, this type of a A- 
term cannot account for more than ~ 0.1 of the critical energy density (see for detailed 



investigation confirming this conclusion). This, unfortunately, prevents us from natural 
explanation of the present A-term with Q\ = (0.5 — 0.8) using "compensation" mechanisms 
or exponential potentials with sufficiently large exponents [T3[ ; in other words, a A-term 



cannot be produced by an exactly "tracker" field as was recently proposed in |TJ . 

A natural and simple description of a variable A-term is just that which was so succes- 
sively used to construct the simplest versions of the inflationary scenario, namely, a scalar 
field with some interaction potential V((p) minimally coupled to the Einstein gravity. Such 
an approach, though phenomenological, is nevertheless more consistent and fundamental 
than a commonly used attempt to decribe a A-term by a barotropic ideal fluid with some 
equation of state. The latter approach cannot be made internally consistent in case of 
negative pressure which is implied by observations [jlO[| , in particular, it generally leads to 
imaginary values of the sound velocity. On the contrary, no such problems arise using the 
scalar field description (this scalar field is called the A- field below). Of course, its effective 
mass \m 2 \ = \d 2 V/d(p 2 \ should be very small to avoid gravitational clustering of this field 
in galaxies, clusters and superclusters. To make a A-term slowly varying, we assume that 
\rrip\ ~ Hq ~ 10 -33 eV, or less (though this condition may be relaxed). Models with a 
time-dependent A-term were introduced more than ten years ago |I5|, and different poten- 



tials V((p) (all inspired by inflationary models) were considered: exponential [T^, 16, 11, 17 
inverse power-law ||18|| , power-law |I9| , cosine pOL Q| . 

However, it is clear that since we know essentially nothing about physics at such energies, 



there exists no preferred theoretical candidate for V(ip). In this case, it is more natural to go 
from observations to theory, and to determine an effective phenomenological potential V(ip) 
from observational data. The two new tests mentioned above are just the most suitable 
for this aim. Really, using the cluster abundance n(z) determined from observations and 
assuming the Gaussian statistics of initial perturbations (the latter follows from the paradigm 
of one-field inflation, and it is in agreement with other observational data), it is possible to 
determine a linear density perturbation in the CDM+baryon dustlike component 5(z) for a 
fixed comoving scale R ~ 8(1 + z)~ x hr l Mpc up to z ~ 1, either using the Press-Schechter 
approximation, or by direct numerical simulations of nonlinear gravitational instability in the 
expanding Universe. 5(z) can be also determined from observation of gravitational clustering 
(in particular, of the galaxy-galaxy correlation function) as a function of z. On the other 
hand, observations of SNe at different z yield the luminosity distance D L (z) through the 
standard astronomical expression m = M+5 log/^+25, where m is the observed magnitude, 
M is the absolute magnitude and D L is measured in Mpc. 

The aim of the present letter is to show how to determine V{<p) from either S(z) or Dl(z), 
and to investigate what additional information is necessary for an unambiguous solution of 
this problem in both cases. The idea has been already annouced by the author in [ft]] , Z2 |, 
now details are given. 

The derivation of V((p) consists of two steps. First, the Hubble parameter H = - = H(z) 
is determined. Here a{t) is the FRW scale factor, 1 + z = a Q /a, the dot means ^ and the 
index denotes the present value of a corresponding quantity (in particular, H(t ) = H(z = 
0) = H ). In the case of SNe, the first step is almost trivial since the textbook expression 
for D L reads: 

D L (z) = a (ri -ri)(l + z), rj = / — - . (1) 



a(t) 



Therefore, 



H(z) = -g- = 'MY 1 
a z arj 




(2) 



Here and below, a prime denotes the derivative with respect to z. Thus, Dl(z) defines H(z) 
uniquely. 

More calculations are required to find H(z) from S(z). The system of background equa- 
tions for the system under consideration is: 

H - — { Pm + T + v ) ' Pm "^G^' (3) 

dV 

(p + 3Hp + — = , 4 

dip 

H = -47rG( Pm + 2 ) . (5) 

Eq. @ is actually the consequence of the other two equations. 

We consider a perturbed FRW background which metric, in the longitudinal gauge (LG), 
has the form: 



ds 2 = (1 + 2$)dt 2 - a 2 {t){l + 2^)5i m dx l dx m } l,m = 1,2,3 



(6) 



The system of equations for scalar perturbations reads (the spatial dependence exp(ikix l ), 
k\k l = k 2 is assumed): 

k 2 

$ = ijr = v , 5 = — -v + 3(v + Hv + Hv) , (7) 
a 1 

$ + H$ = 4:7iG(p m v + ipSif) , (8) 

( k 2 \ dV 

— - + 47rG0 2 $ = 47rG(p m 5 + y5y + 3H<p8<p + — M , (9) 
\ a z ) dip 

( k 2 d 2 V\ dV 
8ip + 3H5ip+[ — + -—)5(p = 4ip$-2—<$>. (10) 

\ or d(p z J dtp 

Eq. (|H]) is the consequence of other ones. Here v and 5<p are, correspondingly, a velocity 
potential of a dustlike matter peculiar velocity and a A-field perturbation in LG, and 5 is a 
comoving fractional matter density perturbation (in this case, it coincides with ( — ) in the 

V P / m 

synchronous gauge). In fact, all these perturbed quantities are gauge- invariant. 

Now let us take a comoving wavelength A = k/a(t) which is much smaller than the 
Hubble radius iJ _1 (t) up to redshifts z ~ 5. This corresponds to A <C 2000 h~ l Mpc at 
present. Then, from Eq. (|10|), 

r a2 /A-^ c dv ^ it- i I dV , a 2 H A 

_( 4v9 $_2 — $), |^| ~ |— <fy>| ~ <Pm|5| . (11) 

Therefore, the A-field is practically unclustered at the scale involved. Now the last of Eqs. (^) 
and Eq. ([|) may be simplified to: 

S=-Kv, -^ = AnGp m 5. (12) 

Combining this with the first of Eqs. (0), we return to a well-known equation for S in the 
absence of the A-field: 

5 + 2H5 - AnGp m 5 = . (13) 

It is not possible to solve this equation analytically for an arbitrary V(ip). Remarkably, 
the inverse dynamical problem, i.e. the determination of H(a) given 6(a), is solvable. After 
changing the argument in Eq. ( |T3"D from t to a (4 = aH-^), we get a first order linear 
differential equation for H 2 (a): 

^^(^t^g-J^. (14) 
da da \ da 2 da J a 6 

The solution is: 

ff2 = 3^ U5\ - * da = + r , z , (15) 

a 6 \da/ Jo da 5' 2 J z 1 + z v 7 

Putting z = in this expression for H, we arrive to the expression of Qq through S(z): 

"° = ^>( 3 f ^tH" ■ (16) 



Of course, observations of gravitational clustering can hardly provide the function S(z) for 
too large z (say, for z > 5). However, 5(z) in the integrands in Eqs. ( |T5|JT6| ) may be well 
approximated by its fl = 1 behaviour (i.e., 5 oc (l + z)" 1 ) already for z > (2 — 3). If massive 
neutrinos are present, one should use here the expression with a written above and with Q\ 
substituted by Q u /Qo (it is assumed that p m includes massive neutrinos, too). 
Finally, using Eq. ([IB]), Eq. ( |T5j ) can be represented in a more convenient form: 

HM_ (i + *>V(o) 3^0!/«i!£U. (17) 



# 2 (0) <5 /2 (;z) u 5' 2 (^) Jo 1 + 2 

Thus, 5(z) uniquely defines the ratio H(z)/H . Of course, appearance of derivatives of S(z) 
in these formulas shows that sufficiently clean data are necessary, but one may expect that 
such data will soon appear. Let us remind also that, for A = const (V((p) = const), we have 

E\z) = Hl{\ - n + no(l + z) 3 ), qo = -1 + = |«o - 1 , (18) 

where go is the acceleration parameter. 

The second step - the derivation of V(ip) from H(a) - is very simple. One has to rewrite 
Eqs. (0,0) in terms of a and take their linear combinations: 

JH 



8nGV(<p) = aH^- + 3H 2 - ^ H 2 f^Y 



and then exclude a from these equations. 

Therefore, the model of a A-term considered in this paper can account for any observed 
forms of Dl(z) and 5(z) which, in turn, can be transformed into a corresponding effective 
potential V(ip) of the A- field. The only condition is that the functions H(z) obtained by 
two these independent ways should coincide within observational errors. Dl{z) uniquely 
determines V(ip), if flo is given additionally (the latter is required at the second step, in 
Eqs. (0)). 5 (z) uniquely determines V(ip) up to the multiplier Hq, the latter has to be 
given additionally to fix an overall amplitude. Observational tests which can falsify this 
model do exist. In particular, a contribution to large-angle CMB temperature anisotropy 
due to the integrated (or, non-local) Sachs- Wolfe effect presents a possibility to distinguish 
the model from more complicated models, e.g., with non-minimal coupling of the A- field to 
gravity or to CDM. However, the latter test is not an easy one, since this contribution is 
rather small and partially masked by cosmic variance. 
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